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Summary. Aldous and Pitman (1994) studied asymptotic distributions as n — > oo, 
of various functionals of a uniform random mapping of the set {1, . . . , n}, by con- 
structing a mapping-walk and showing these random walks converge weakly to a 
reflecting Brownian bridge. Two different ways to encode a mapping as a walk lead 
to two different decompositions of the Brownian bridge, each defined by cutting the 
path of the bridge at an increasing sequence of recursively defined random times in 
the zero set of the bridge. The random mapping asymptotics entail some remark- 
able identities involving the random occupation measures of the bridge fragments 
defined by these decompositions. We derive various extensions of these identities 
for Brownian and Bessel bridges, and characterize the distributions of various path 
fragments involved, using the Levy-Ito theory of Poisson processes of excursions for 
a self-similar Markov process whose zero set is the range of a stable subordinator of 
index a € (0, 1). 
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1 Introduction 

In a previous paper pP we showed how features of a uniformly distributed 
random mapping M„, from [n] := {1, 2, . . . , n} to itself, could be encoded as 
functionals of a particular non-Markovian random walk on the non-negative 
integers. This mapping-walk, suitably rescaled, converges weakly in C[0, 1] as 
n — > oo to the distribution of the reflecting Brownian bridge defined by the ab- 
solute value of a standard Brownian bridge B^'^ with Bq'^ = B^'^ = obtained 
by conditioning a standard Brownian motion B on Bi = 0. Two important 
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features of a mapping are the vector of sizes of connected components of its di- 
graph, and the vector of sizes of cycles in its digraph. Results of imply that 
for a uniform random mapping, as n — > oo, the component sizes rescaled by n, 
jointly with corresponding cycle sizes rescaled by y^, converge in distribution 
to a limiting bivariate sequence of random variables (A/^ , )j=i.2.... where 
(/j)j=i,2,... is a random interval partition of [0,1], with A/, the length of Ij 
and Lj . the increment of local time of B''" at over the interval Ij . With the 
convention for ordering connected components of the mapping digraph used 
in T, the limiting interval partition is {Ij) = {if): according to the following 
definition. Here, and throughout the paper, t/, [/i, C/2, • ■ ■ denotes a sequence 
of independent uniform (0, 1) variables, independent of B^^ . 

Definition 1 (the D-partition [J). Let Ij^ :— [Dvj_i, Dvj] where Vo = 
Dyg — and Vj is defined inductively along with the Dy^ for j ^ 1 as 
follows: given that Dv^ and Vi have been defined for ^ i < j, let 

Vy.^Dv^_,+U,{l^Dv,_J, 

so Vj is uniform on [Dy^_-^,1] given B^'^ and {VijDy.) for ^ i < j, and let 

Dv, :=inf{t ^ V.J : Sf"' = 0}. 

On the other hand, a variation of the main result of 1 shows that with 
a different ordering convention, the mapping component sizes rescaled by n, 
jointly with their cycle sizes rescaled by ^/n, have a limit distribution speci- 
fied by the sequence of lengths and Brownian local times (A/^ , )j=i,2.... a 
differently defined limiting interval partition. This is the partition {Ij) — {ij) 
defined as follows using the local time process (L" , ^ w ^ 1) of B'°'^ at 0: 

Definition 2 (the T-partition). Let ij := [Tj-i,Tj\ where Tq := 0, Vo := 
0, and for j ^ 1 

% ■■= 1 - - U^). (1) 

1=1 

so Vj is uniform on [Vj-i, 1] given B^" and {Vi,Ti) for ^i < j, and 
Ty.= inf{u:Ll/L°,>V,}. 

For each of these two random interval partitions {Ij) we are interested 
in the distribution of the bivariate sequence of lengths and local times 
{Xj. , ^7^ )^=1, 2,... and the distribution of the associated path fragments B^'^[Ij] 
and standardized fragments i3^''[/j]. Here for a process X := {Xt,t e J) pa- 
rameterized by an interval J, and / = [G/, Dj] a subinterval of J with length 
A/ := Dj — G/ > 0, we denote by X[I] or X[Gi, Dj] the fragment of X on I, 
that is the process 
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:= Xg,+u (0 < « < A,). (2) 

We denote by X*[7] or X*[Gi,Di] the standardized fragment of X on 7, 
defined by the Brownian scaling operation 

X[I]u\i Xqi+uXi 



/AT 



(0<u< 1). 



(3) 



Figure 1 illustrates these definitions for a typical path oi X = B^^. Note 
that the first interval if of the -D-partition ends at the time of the 
first zero of B^'^ after a uniform(0, l)-distributcd time Ui, whereas the first 
interval Ii of the T-partition ends at the time Ti when the local time of 
B^'^ at has reached a uniform(0, l)-distributed fraction of its ultimate value. 
As illustrated in Figure 1, the associated fragments of B^'^ arc qualitatively 
different: B^^[iP] ends with an excursion while B^'[ll] does not. 
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Figure 1. The two interval partitions. 
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Despite this difference between tfie fragments of B over the D- and T- 
partitions, the random mapping asymptotics have the following corollary. Let 
(/^^j) and denote the length-ranked D-partition and the length-ranked 

T-partition respectively, meaning /^^^ is the jth longest interval in the D- 
partition, and l'^^^ is the jth longest interval in the T-partition. 

Theorems. Considering the four bivariate sequences (Xi^ , L^.)j=i^2.... of 
lengths and bridge local times at 0, for (Ij) one of the four random inter- 
val partitions of [0,1] defined by {Ij'), (-^(^)); i^J) '^'^ (^(j))' 

(i) the bivariate sequence has the same distribution for (^(^-j) as for {ifj-^)', 

(ii) the bivariate sequence for {if') is the bivariate sequence for {I^j^) in a 
length-biased order; 

(iii) the bivariate sequence for {ij) is the bivariate sequence for {I'^j^) in an 
L^-biased order; 

(iv) the sequence of local times {L^.) has the same distribution for {Ij') as 
for {ij), whereas the sequence of lengths {\i-) does not. 

See ESI for background about size-biased random orderings. To illus- 
trate the meaning of (iii) for instance, for each k ^ 1, conditionally given 
the entire bivariate sequence {Xjt , L^,j. )j=i 2 , the probability of the event 

{I^ — ^l^)) is L^jT /Li, where ~ J2j I^'jt almost surely. And given 

also {I-f = ^(^-j), for each m ^ 1 with m ^ k the probability of the event 

(/J — ijly^-^) is / (^1 ^ -^/r) ^^"^ ^'^ another way, parts (i)-(iii) of 

the corollary state that the bivariate sequence {\i. , L^j. )j=i,2,... for {Ij) — {Ij^) 
is distributed like a length-biased rearrangement of the bivariate sequence for 
{Ij) = {iJ"), which is in turn distributed like an i^-biased rearrangement of 
the bivariate sequence for {Ij) ~ {if)- Consequently, the distribution of any 
one of the four bivariate sequences determines the distribution of each of the 
others. 

The rest of this paper is organized as follows. Section |21 explains how we 
discovered Theorem |21 by consideration of random mapping asymptotics. We 
recall the theorem from 1 which describes the asymptotics of mapping-walks 
in terms of the fragments of B^'' defined by the Z?-partition, and present the 
companion result, for a different ordering of components, where the limit in- 
volves the fragments of B^'^ defined by the T-partition. Section O lays out 
our results regarding the decomposition of B^'' into path fragments associ- 
ated with the D- and T-partitions, in a way which does not depend on the 
random mapping asymptotics. In particular, we describe the three different 
distributions of bivariate sequences featuring in the three parts of Theorem^l 
We formulate and prove these results more generally, for B^^ the standardized 
bridge of a recurrent self-similar Markov process B whose inverse local time 
process at is a stable subordinator of index a for some a € (0, 1). So a = 1/2 
for B a standard Brownian motion as supposed in previous paragraphs, and 
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a — 1~ 5/2 for B a Bessel process of dimension S £ (0, 2). Some of the results 
in Section|31 like Theorem|3| can be viewed in the Brownian case as asymptotic 
counterparts (under weak convergence of mapping- walks) of some combinato- 
rial symmetries of random mappings, discussed in Section f2.2l Other results in 
the Brownian case, especially those involving the method of Poissonization by 
random scaling 1341 OHl , are not obvious from the combinatorial perspective, 
but provide explicit limit distributions for functionals of uniform random map- 
pings. See also [Sj where we apply this method to characterize the asymptotic 
distribution of the diameter of the digraph of a uniform mapping. Sections 0] 
and|5lprovide some proofs and further details of the main results in Section^ 
while Section |5| contains various complements. In particular, we show in Sec- 
tion l6.2l that Theorem 13 holds even more generally for interval partitions {Ij^) 
and (ij) defined as before, but with the random zero set of B^' replaced by 
the complement of IJ^- Ij^, where 1^^ is any exchangeable random partition of 
[0, 1] into an infinite number of intervals, and (L° , ^ u ^ 1) is the associated 
local time process, as defined by Kallenberg jl6j. This is the limiting case of a 
corresponding result for a finite exchangeable interval partition of [0, 1], which 
we prove by a combinatorial argument. 

In companion papers |2] and |S] we show that Brownian bridge asymptotics 
apply for models of random mappings more general than the uniform model, in 
particular for the p-mapping model j24ll29| . and that proofs can be simplified 
by use of Joyal's bijection between mappings and trees. See also (30) for a 
recent review of the applications of Brownian motion and Poisson processes 
to the asymptotics of various kinds of large combinatorial objects, including 
partitions, trees, graphs, permutations, and mappings. 

2 Random Mappings 

In this section we explain how study of random mappings led us to consid- 
eration of the two interval partitions of Brownian bridge, and show how the 
distributions of path fragments of the bridge defined by these partitions en- 
code various asymptotic distributions for mappings. 

2.1 Mapping- walks and the two orderings 

A mapping M„ : [n] [n] can be identified with its digraph of edges 
{(i, Mn{i)), i G [n]}. The connection between random mappings and Brown- 
ian bridge developed in T can be summarized as follows. 

• A mapping digraph can be decomposed as a collection of rooted trees to- 
gether with extra structure (cycles, basins of attraction). 

• A rooted tree can be coded as a discrete tree-walk, a walk excursion starting 
and ending at 0. 
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• Given some ordering of tree-components, one can concatenate walk-excur- 
sions to define a discrete mapping-walk which codes M„. 

• For a uniform random mapping, the induced distribution on tree-compo- 
nents is such that the tree- walks, suitably normalized, converge to Brownian 
excursion as the tree size increases to infinity. 

• So for a uniform random mapping, we expect the mapping-walks, suitably 
normalized, to converge to a limit process defined by some concatenation 
of Brownian excursions. 

• With appropriate choice of ordering, the limit process is in fact reflecting 
Brownian bridge. 

We now amplify this summary, emphasizing the only subtle issue - the 
choice of ordering. Fix a mapping M„. It has a set of cyclic points 

Cn ■■= {i e [n] : M^{i) = i for some k ^ l}, 

where is the fcth iterate of Af„. Let Tn^c be the set of vertices of the 
(perhaps trivial) tree component of the digraph with root c G C„. The tree 
components are bundled by the disjoint cycles Cn,j Q Cn to form the basins 
of attraction (connected components) of the mapping digraph, say 

^n,j [J %i,c^C„j with US„j = [n] and yC„j=C„ (4) 

cGC„,j j j 

where all three unions are disjoint unions, and the Bnj and Cnj are indexed in 
some way by j = 1, . . . , Kn say. The construction in encodes the restriction 
of the digraph of M„ to each tree component Tn^c of size k (that is, with k 
vertices) by 2k steps of a tree-walk with increments ±1 on the non-negative 
integers. The tree-walk proceeds by a suitable search of the set Tn,c, making 
an excursion which starts at and returns to for the first time after 2k 
steps, after reaching a maximum level 1 -I- hn{c), where hn{c) is the maximal 
height above c of all vertices of the tree Tn^c with root c, that is 

K{c) = max{/i : 3i(^[n] with M^{i) = c and Ml{i) ^ C„ for j < /i}. (5) 

It was shown in j4j that as fc — > cx), the distribution of the tree- walk for a 
/c-vertex random tree, of the kind contained in the digraph of the uniform 
random mapping Af„ for k ^ n, when scaled to have 2k steps of ±l/^/k per 
unit time, converges to the distribution 2B°^ for a standard Brownian 
excursion. Subsequent work P2^ shows that the same result holds for a variety 
of codings of trees as walks. Consequently, any of these codings would serve 
our purpose in the following definitions. 

We now define a mapping- walk (to code M„) as a concatenation of its tree- 
walks, to make a walk of 2n steps starting and ending at with exactly |C„| 
returns to 0, one for each tree component of the mapping digraph. To retain 
useful information about Af„ in the mapping-walk, we want the definition of 
the walk to respect the cycle and basin structure of the mapping. Here are 
two orderings that do so. 
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Definition 4 (cycles-first ordering). Fix a mapping Af„ from [n] to [n]. 
If Mn has Kn cycles, first put the cycles in increasing order of their least 
elements, say c„^i < c„ 2 < • • ■ < Cn^K^- Let Cn,j be the cycle containing Cnj, 
and let Bn,j be the basin containing Cn,j ■ Within cycles, list the trees around 
the cycles, as follows. If the action of Mn takes Cnj —>•••—!■ Cn.j 

for each 1 ^ j ^ Kn, the tree components Tn^c o-fs listed with c in the order 

(Cn,ia, . . . , Cn,l, C„,2,l, • ■ • , Cn,2, ■ ■ ■ , Cn,K„,l, ■ ■ ■ , Cn, (6) 

The cycles-first mapping-walk is obtained by concatenating the tree walks de- 
rived from AIn in this order. The cycles-first search of [n] is the permutation 
a : [n] [n] where aj is the jth vertex of the digraph of Mn which is visited 
in the corresponding concatenation of tree searches. 

Definition 5 (basins-first ordering [1 ). // M„ has Kn cycles, first put 
the basins Bn,j in increasing order of their least elements, say 1 = < 
bn.2 < . . . < bn,K„; let Cnj G Cn,j be the cyclic point at the root of the tree 
component containing bnj- Now list the trees around the cycles, just as in 
but for the newly defined Cnj and Cnj^i- Call the corresponding mapping-walk 
and search of [n] the basins-first mapping-walk and basins-first search. 

Be aware that the meaning of Bnj and Cnj now depends on the ordering 
convention. Rather than introduce two separate notations for the two order- 
ings, we use the same notation for both, and indicate nearby which ordering is 
meant. Whichever ordering, the definitions of Bnj and Cnj are always linked 
by Bnj 2 C„j , and holds. 

Let us briefly observe some similarities between the two mapping- walks. 
For each given basin B of M„ with say b elements, the restriction of M„ to 
B is encoded in a segment of each walk which equals at at some time, and 
returns again to after 2b more steps. If the basin contains exactly c cyclic 
points, this walk segment of 26 steps will be a concatenation of c excursions 
away from 0. Exactly where this segment of 26 steps appears in the mapping- 
walk depends on the ordering convention, as does the ordering of excursions 
away from within the segment of 26 steps. However, many features of the 
action of M„ on the basin B are encoded in the same way in the two different 
stretches of length 26 in the two walks, despite the permutation of excursions. 
One example is the number of elements in the basin whose height above the 
cycles is h, which is encoded in either walk as the number of upcrossings from 
/i to /i + 1 in the stretch of walk of length 26 corresponding to that basin. 

2.2 Symmetry properties of random mappings 

We now apply the definitions above to a uniform random mapping M„. Of 
course, the random partition {/Bn j}j=i,...,i<'„ of [n], and the random partition 
{Cn.j}j=i,...,Kn of Cm are the same no matter which ordering convention is 
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used. Each random partition is exchangeable, meaning its distribution is in- 
variant under the action of a permutation of [n] . Let us spell out some further 
symmetry properties, each of which turns out to have some analog in the 
limiting Brownian scheme. 

(a) The cycles-first ordering has the following very strong symmetry prop- 
erty: conditionally given \Cn\ = m the tree components in cycles- first ordering 
form an exchangeable sequence of m random subsets of [n] ; moreover this ex- 
changeable sequence is independent of the sequence of cycle sizes |C„j | with 
X]j l^n.il = Consequently, given |C„| — m, the cycles-first mapping- walk is 
a concatenation of m exchangeable excursions away from 0, and this mapping- 
walk is independent of |C„j|, j — 1,2, ... , Kn. 

(b) The basins-first ordering does not share the symmetry property above. 
But it has a different one: given that the basin Bn,i containing 1 has size 
\Bn,i\ = b, the action of M„ on [n]—Bn,i is that of a uniform random mapping 
of a set oi n — b elements. So given \Bn,i\ = b, the basins-first mapping- walk 
decomposes after 2b steps into two independent segments: the first 2b steps are 
distributed like the basins-first walk for a uniform mapping of [b] conditioned 
to have a single basin, and the remaining 2{n — b) steps distributed like the 
basins-first walk associated with a uniform mapping of [n — b] . 

(c) The sequence of basin sizes 1 ^ :/ ^ Kn) does not have the same 
distribution for both orderings. For instance, if |;B„^i| = 1 in the basins-first 
ordering then \Bn,i\ = 1 in the cycles-first ordering, but (for n ^ 3) not 
conversely. So the distribution of \Bn,i \ must be different in the two orderings. 

(d) For a given mapping M„, the sequence of cycle sizes (|C„j |, 1 ^ j ^ Kn) 
may be different for the two different orderings. But for M„ with uniform 
distribution on the two sequences of cycle sizes have the same distri- 
bution: given \Cn\ = m, either sequence is distributed like the sizes of cycles 
of a uniform random permutation of [m] in the (size-biased) order of least 
elements of the cycles. That is to say, given |C„| — m, the distribution of 
\Cn,i\ is uniform on [to]; given |C„| — to and \Cn,i\ with |C„| — \Cn,i\ = mi, 
the distribution of \Cn,2\ is uniform on [toi], and so on. This is a well known 
property of uniform random permutations for the cycles-first ordering, and 
was shown for the basins-first ordering in fP, Lemma 22]. 

2.3 Brownian asymptotics for the mapping-walks 

We now come to the main point of Section El the definitions of the interval 
partitions of Brownian bridge are motivated by the following theorem. 

Theorem 6. The scaled mapping-walk (Afi"',0 ^ u ^ 1), with 2n steps of 
±1/y^ per unit time, for either the cycles-first or the basins-first ordering 
of excursions corresponding to tree components, converges in distribution to 
2\B'°''\ jointly with 
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where (i" , ^ m ^ 1) is the process of local time at of B^'^ , normalized so 
that P{Ll > l)^ e~^'/2. Moreover, 

(i) for the cycles-first ordering, with the cycles Bn,j in order of their least 
elements, these two limits in distribution hold jointly with 

as j varies, where the limits are the lengths and increments of local time of 
B^^ at associated with the interval partition (Ij) whereas 

(ii) Q basins-first ordering, with the basins Bn,j listed in order of 
their least elements, the same is true, provided the limiting interval partition 
is defined instead by (Ij) := (-^j')- 

The result for basins- first ordering is part of ^ Theorem 8]. The variant 
for cycles-first ordering can be established by a variation of the argument in 
PP, exploiting the exchangeability property of the cycles-first ordering (Sec- 
tion |^fa) instead of Section fbll. See also JO] and [S] for alternate ap- 
proaches to the basic result of p. 

We now explain how we first discovered some of the facts about Brown- 
ian bridge presented in Theorem O by consideration of Theorem |21 and the 
symmetry properties of Section [2. 21 The arguments below are not part of the 
formal development in this paper. Indeed we show in Section that the re- 
sults of Theorem hold much more generally, so these results do not really 
involve much of the rich combinatorial structure of mapping digraphs involved 
in Theorem El 

(a) In the basins-first ordering, the first basin is by definition the basin con- 
taining element 1, and its walk-segment ends at the first time that the walk 
returns to after the basins-first search has reached element 1. Suppose we 
could replace element 1 by a uniform random element, so the walk-segment 
corresponds asymptotically to the walk-segment ending at the first time of 
reaching after a uniform random time on [0, 2n\. Rescaling, this corresponds 
to the time interval [0, DyJ in Definition Of course it is not obvious, and 
indeed is somewhat counter-intuitive, that replacing element 1 by a uniform 
random element will preserve length of walk-segment. But this is what even- 
tually will emerge from our calculations. 

Now consider the cycles-first ordering. The first basin is by definition the 
basin containing the smallest-numbered cyclic element c„.i, and its walk- 
segment ends at the first time after reaching element c„.i that the walk re- 
turns to 0. Suppose as before (and again this is not obvious) one can replace 
element Cn,i by a uniform random cyclic element, so the walk-segment corre- 
sponds asymptotically to the walk-segment ending at the first time of reaching 
after visiting C/*|C„| cyclic vertices, where U* has uniform[0, 1] distribution. 
Rescaling, this corresponds to the time interval [0,ri] in Definition[21 
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(b) The recursive property of the basins-first ordering in Section r2.2f b'l plainly 
corresponds, under the asymptotics of Theorem El to the recursive decompo- 
sition of Brownian bridge at time Dy-^ described later in Lemma |51 

(c) In Section I^Jc) we observed that the distribution of Bn.i was different 
in the two orderings. This difference persists in the limit: Theorem |B1 and the 
calculation below H26|) imply 



J E{Dvi) = 2/3 (for the basins-first ordering) 
[ E{Ti) — 1/2 (for the cycles-first ordering). 



(d) It is well known 01] that the asymptotic distribution as n ^ oo of the 
fractions of elements in cycles of a random permutation of [n] , with the cycles 
in order of their least elements, (which amounts to a size-biased random order 
by exchangeability), is the uniform stick-breaking sequence Uj Y[i=i (1 ~ Ui). 
So the convergence in distribution Q of \Cn\l\pri to L\, and the "uniform 
random permutation" feature of the cyclic decomposition (Section |^2Id)), 
combine to show that with either ordering \Cn.j\l \fn V\. with the same 
joint distribution: 

^ 1) - (^?^/nV - v.). 3 ^ (9) 

for both Ij — and Ij = ij. This is part (iv) of Theorem O which is 
generalized later by l|77|l and Theorem 

(e) Let Bn^(j) be the jth largest basin of M„, with some arbitrary convention 
for breaking ties, and let C„ (j) be the cycle contained in Bn.(j)- It follows 
immediately from the convergence in distribution ((SJ that 



\Bn,{j)\ \Cn,{j)\ 



(10) 



jointly as j varies, where /(j) is the length-ranked interval partition derived 
from either {if) or [ij). This is part (i) of Theorem O By exchangeability 
considerations, before passage to the limit the bivariate sequence in (jSJ as j 
varies is that in H10|l biased by cycle-size in the cycles-first order and biased 
by basin-size in the basins-first ordering. Hence the conclusions of parts (ii) 
and (iii) of Theorem |21 by a straightforward passage to the limit, 
(f) Due to Section ITW a). it makes no difference to anything if in the cycles- 
first ordering we replace the ordering within the jth cycle Cn,j,i, Cn,j,2, • • ■ , c„j 
by the possibly more natural c„j-, c„j^i, c„j-_2, • ■ ■ , c„ j But in the 

basins-first ordering, this innocent looking change would spoil convergence 
to 2|_B'"'|. This is because in the basins-first ordering the tree with root c„_i 
is the tree containing 1, which is a size-biased choice from the exchangeable 
random partition of [n] into tree components. As such, it tends to be a big 
tree. In fact, results from ^ imply that, if the mapping- walk is started by 
the excursion coding the tree rooted at c„^i, the limit process will start with 
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a zero free interval whose length is distributed as Djj — Gu in Lemma El be- 
low for a — 1/2. Such a process is obviously not 2|i?'^''| or any other familiar 
Brownian process. 

(g) The proof of Theorem yields more information about the asymptotic 
sizes of tree components than can be deduced from the statement of that 
theorem. For instance, if are the ranked sizes of the tree components 

of Mn, and Hn,i are the corresponding maximal tree heights, as in 0, then 
HuaI \/n)i=i_2,--- converges in distribution to the sequence of ranked 
lengths and corresponding maximal heights of excursions of 2|_B'^''|, whose 
distribution was described in 35, Theorem 1 and Example 8]. If only the 
tree components of Bn,j were considered, the limit would be derived from 
excursions of B'"^ over the appropriate random interval Ij as in Theorem |H1 
with joint convergence as j varies. 

3 The bridge decompositions 

This section presents our main results for the D- and T-partitions. For ease 
of comparison, the results are presented together here, with outlines of the 
proofs. Some proofs and further details are deferred to Section 0] for the D- 
partition, and to Section [S] for the T-partition. Our primary interest is the 
analysis of the D- and T-partitions derived from a standard Brownian bridge, 
and the connections between these random partitions and the asymptotics of 
random mappings discussed in Section|21 But we find that our analysis applies 
just as well to the D- and T-partitions for a standardized bridge B^^ derived 
from B a recurrent self-similar Markov process whose inverse local time pro- 
cess at is a stable subordinator of index a for some a S (0, 1). Readers who 
don't care about this generalization can assume throughout this section that 
B is standard one-dimensional Brownian motion, and a = (3 ^ 1/2. 

3.1 General framework 

Following Pitman- Yor 1351 §2], we make the following basic assumptions: 

• B :— {Bt,t ^ 0) is a real or vector- valued strong Markov process, started 
at Bo — 0, with state space a cone contained in M.'^ for some d = 1,2,..., 
and cadlag paths. 

• -B is /3-self-similar for some real f3. That is to say, if i?, [0, t] now denotes 
the standardized process derived from B on [0, t] as in (0, using \j instead 
of ^/Xi in the denominator, then B^[0,t] = B[0, 1] for all t > 0. 

• The point is a regular recurrent point for B, meaning that almost surely 
both and oo are points of accumulation of the zero set of B. 

As a well known consequence of these assumptions E^l, there exists a 
continuous local time process for B at 0, say {L^{B),t ^ 0), whose inverse 
process 
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Ti ~ inf{t : ^{B) > e\ 0) 
is a stable subordinator of index a for some a E (0, 1). That is 

EeM-^ri)^expi-icC) > 0) (H) 

for some c > 0, in which case 

LOJB) = '^^^'"^ lim e^Nt (12) 

uniformly for bounded t almost surely, where Nt g{B) is the number of excur- 
sion intervals of B in [0, t] whose length is greater than e. Formula (|12|l can 
be then used with X instead of B to define L^X) for various other processes 
X derived from B by conditioning or scaling, such as the standardized bridge 
B^'^ introduced in the next paragraph. As a consequence of H12() with X in- 
stead of i?, there is following basic a-scaling rule for such local time processes: 
for / = [Gi, Di] a random subinterval of length A/ := Di — Gi contained in 
the time domain of X, and L%X) := L%^{X) - L%^{X), 

L%X) = \JL\{X^[I]). (13) 

Associated with the self-similar Markov process B are corresponding distribu- 
tions of a standard B-bridge B^^ , a standard B-excursion i?™, and a standard 
B -meander 5™°, defined by the following identities in distribution, valid for 
all t > 0: 

Gt] = S'^'-; B4Gt,Dt]=B'''; [G* , = (14) 

where Gf :— Gt{B), Dt := Dt{B), and for any process X we use the notation 

Gt{X) -.^ sup{u < t : Xu ^ 0} (15) 
DtiX) := inf{u ^t: X^ = 0}. (16) 

See for a review of properties of B^'^, B°^ and B™'^ in the Brownian case 
when B is Brownian motion with state space M, and /? = a = 1/2. See §3] 
and |34j for some treatment of i?*^' and B'^^ in the Bessel case when B with 
state space M^o is a recurrent Bessel process of dimension (5 = 2 — 2a G (0, 2), 
and /3 = 1/2. Other examples are provided by recurrent stable Levy processes 

symmetrized or skew Bessel processes and Walsh processes |S1[7|. 

According to the Levy-Ito theory of excursions of B, applied to the stan- 
dard B-bridge as in if (Ij^) is the interval partition of [0, 1] defined 
by the excursion intervals of B^'^ in length-ranked order, then the processes 
B^'^[Ij^] are i.i.d. copies of B'^^, independent of (Ij^), which is an exchange- 
able interval partition in the sense of |16) recalled in Section 16.21 Moreover, 
the distribution of ranked lengths (A/™) depends only on a, as described in 
123 (16)] and Example 8]. This general excursion decomposition of B^'^ 
implies that various results known for Bessel bridges hold also in the present 
general setting, and we take this for granted without further comment. 
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3.2 Main Results 

All results of this section are presented with the notation and general frame- 
work of the previous section: B^'^ is the standard B-bridge derived from a 
self-similar recurrent Markov process B whose continuous local time process 
{L^{B),t ^ 0) is the inverse of a stable subordinator {ti,£ ^ 0) of index 
a e (0, 1). The D- and T-partitions are defined in terms of B^' and its local 
time process at 0, according to Definitions ^ and 13 

Theorem |3| presented in the introduction in the Brownian case, is true in 
the more general framework of this section, as a consequence of the following 
theorem: 

Theorem 7. For a random interval I C [0,1], let fij denote the random 
occupation measure induced by the path of B^^ll], so for each Borel subset A 
of the state space of B^'^ 

^iIiA) := jl{Bf eA)dt. 

(i) The sequence of occupation measures {fJ.i ) has the same distribution for 
each of the two length-ranked partitions [Ij ) — (/^^^ ) and [Ij ) = i-^fj) ) ■ 

(ii) For {Ij) = {if) the sequence of occupation measures {^i-) is in X-biased 
order, where A/ . is the total mass of the random measure fij, . 

(iii) For {Ij) = {ij) the sequence of occupation measures (/i/^) is in -biased 
order, where L^, := L^j.{B^''). 

(iv) For {Ij) equal to any one of the four interval partitions {if), i-^lj)), i^J) 
or {ijj-^), conditionally given Xj. — Xj and L^, — ij for all j ~ 1,2,..., 
the random occupation measures ^i., j = 1,2,... are independent, with ^j. 
distributed like the random occupation measure of a process with the common 
conditional distribution of 

(B[0, t]\Bt= 0, = £) = (B[0, n] \n^t) (17) 
for t — Xj and I — ij. 

Proof. Propositions llUllTTI andiniprovide more explicit descriptions of the law 
of (A/j , Lj. , i3j?''[/j])j=i^2,..., for each of the four interval partitions {Ij). The 
above results for occupation measures are deduced from these propositions 
using Lemma El The fundamental switching identity \V7\ is well known \M\ 
§5]. □ 

By general theory of local time processes for diffusions or continuous semi- 
martingales ^1 EHl ISH] ! in the Brownian and Bessel cases for each random 
subinterval / of [0, 1] the random occupation measure /x/ derived from B'"^ 
has an almost surely continuous density relative to m at cc, where m is a 
multiple of the speed measure of the one-dimensional diffusion B. To be precise 
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about normalization of local times, in the Brownian case with state space M, 
we take m{dx) — dx, so that Hll|l holds with a — 1/2 and c = \/2. In the 
Bessel((5) case with state space M^o, we take m{dx) = 2x^~^ dx, so that Hll|l 
holds with a = 1 - (5/2 and c = 2}-"r{l - a)/r(a), by |S1 (7.c)]. In either 
case, Lj, in (iii) and (iv) is recovered like A/^. as a measurable function of the 
random occupation measure fij^ . The distribution of the local time density of 
the conditional occupation measure in (iv) is described by a conditional form 
of the Ray-Knight theorem: see |2()l |2H| for details in the Brownian case. 

Our analysis of the D-partition is the following expression of the decom- 
position of B^'^ at the times Dy^ , implicit in P in the Brownian case: 

Lemma 8 (U). For each j, the pre-Dvj fragment of the bridge i?^' [0, Dy,] 
is independent of the standardized post-Dy fragment B^^ [Dv ,1], which has 
the same distribution as B^' . 

This is easily verified, because the Dy^ are stopping times relative to a 
filtration with respect to which B^'^ has a strong Markov property. 

To describe various distributions, let (Fg, s ^ 0) denote a standard gamma 
process, that is the increasing Levy process with marginal densities 

P{F, £ dx)/dx = x'~^ e"^ {x > 0), (18) 
1 [sj 

so Ft — Fg = Ft-s for < s < t. Recall that for a,b > the beta(a, h) 
distribution is that of 

Pa,b '■— Fa/Fa+b, which is independent of Fa+b, with (19) 



PWa.b e du) = ) l> u^-\l ~ uf-^ du (0<«<1). (20) 
F{a)F{b) 

It is well known [251 Lemma 3.7] that for Gt = Gt{B), 

the standard i?-bridge B^[0, Gt] is independent of Gt with Gt/t = i~a- 

' (21) 

Lemma 9 (ID Prop. 2], |23 Prop. 15]). Let U with uniform[0,l] dis- 
tribution be independent of B^'^ , and let Gjj := Gu{B^'^), Djj := Dij{B^'^). 
Then 

{Gu,Du — Gu, 1 — F>u) = {Fa, A — Fa,Fi+a — A)/A+a- 

Moreover, the random vector {Gu, Du — Gu, 1 — Du) and the three standard- 
ized processes B^'[0,Gu], B^'[Gu,Du] and B^'^lDu,!] are independent, with 



B^'[0, Gu] = B^'[Du, 1] = B^' and B^'[Gu, Du] = B'^''. (22) 
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Proposition 10. For the D-partition 

(i) the sequence of lengths is such that 

Xjn ^WjY[{1~W^) (23) 

1=1 

for a sequence of independent random variables Wj with Wj = 

(ii) The corresponding sequence of local times at can be expressed as 

L%=X'}.L%B^^[I^]) (24) 

3 J 

where the Li{B^'^[Ij^]) are independent random variables, independent also of 
the lengths Xjp , with 

L?(B^'-[/f])^L°(B)^rr" (25) 

for Ti with the stable distribution of index a defined by (|11|) . 

(iii) The standardized path fragments B^'^[Ij'] are independent and identically 
distributed like -B^''[0, and independent of the sequence of lengths (Xjd). 

(iv) For the length-ranked D-intervals I^^ instead of Ij^ , the lengths (Xio^) 
have the Poisson-Dirichlet{a) distribution defined by ranking (A/u) as in (i), 
while parts (ii) and (iii) hold without change. 

Proof. Parts (i)-(iii) are obtained by repeated application of Lemmas |H1 and 
using the a-scaling rule H13(l for local times and l|21(l . as in |25l Lemma 3.11], 
for part (ii). The second identity in (ii) is a well-known consequence of the 
inverse relation between {L^{B),t > 0) and (r^,^ ^ 0), as discussed in [ST]. 
Part (iv) follows immediately from (i)-(iii). □ 

See ^1 021 and Lemma El for background on the Poisson-Dirichlet dis- 
tribution appearing in (iv). Levy |21| showed that in the Brownian case the 
common distribution of L\{B) and appearing in (|25|l is simply the dis- 

tribution of with Bi standard Gaussian. But this does not generalize to 

the Bessel((5) case for general 5 = 2 — 2a. Then Bi = ^J2F\Z^^ which is a 
simple transformation of the stablc(a) distribution of t\ only for a = 1/2. 

The difficulty involved in Theorem is that Definition (21 of the Tj involves 
the local time L\ := L]'(i?'^'), which depends on the path of B^"^ over the 
whole interval [0,1]. While we can describe the finite-dimensional distribu- 
tions of the bivariate sequence (A/t, L°r)j=i,2,... by conditioning on L\ (see 

Proposition I20|l . this description is more complicated than our description of 
{XjD , L^,D)j=i,2.... in Proposition [TT)I In particular, 

XjT^Xjn. (26) 
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Indeed, by we have 

E{Xjp) = E{Wi) = 1/(1 + a) > 1/2, 

whereas (by symmetry of B^'^ with respect to time reversal in the Brownian 
or Bessel case) the distribution of Ti is symmetric about 1/2, so whatever 
a e (0,1) 

E{Xjt) ^ EiTi) - 1/2. 

Still, as explained combinatorially in the Brownian case around lO, the two 
partitions give rise to the same distribution for the sequence of local times: 




where the second equality by definition is read from (^). The first equality 
in distribution of sequences follows from Lemma [S] and the consequence of 
Lemma|21 noted in [Jj (3)-(4)] in the Brownian case, that 

L'^jo/LI has uniform distribution on (0, 1), and is independent of L^. (28) 

As indicated in Section 16.21 this can also be checked in general using the 
exchangeability of the excursion interval partition. 

According to Proposition^! the standardized bridge fragments over in- 
tervals of the Z?-partition are i.i.d. copies of B^^lO, Djj], both for the intervals 
in their original order and for the intervals in length-ranked order. A subtle 
feature of the T-partition is that the standardized bridge fragments over its 
intervals are neither independent nor identically distributed in their original 
order, but these fragments become i.i.d. when put into length-ranked order. 
This and other parallels between the T- and D-partitions in length-ranked 
order are presented in the following Proposition: 

Proposition 11. For the T-partition in length-ranked order 

(i) the sequence of lengths {Xjt ) has the same Poisson-Dirichlet{a) distri- 

hution as {Xjd ). 

(3) 

(ii) The corresponding sequence of local times at can be expressed as 

L?. =A?. i;(S-[/5,]) (29) 

where the L5(-^*'^[-^(j)]) '^'"^ independent random variables, independent also 

of the lengths L'It , with 

if) 

Lm'[lu)])=L'i{B)^rr (30) 

just as in (|25|) . 

(iii) The standardized path fragments B^'^ [/|^^ ] are independent and identically 
distributed like _B,[0,ti], and independent of the sequence of lengths (Ajr ). 
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The only difference between this description of the law of the sequence 
{Xij , L'j B^'^[Ij])j^i 2.... for Ij — ij-.y and the previous description in Propo- 
sition^] for Ij — I-^^y is that the common distribution of the standardized T- 
fragments is that of B^^[Q^ Djj]^ whereas the common distribution of the stan- 
dardized D-fragments is that of B*[0,ri]. The standardized process i?*[0,ri] 
is known as the pseudo-bridge associated with the self-similar Markov process 
B. The following Lemma was established by Biane, Le Gall and Yor ^I] in 
the Brownian case, and extended to the Bessel case in [21 Theorem 5.3]. 

Lemma 12 ([ 111 EH)- The law of the pseudo-bridge _B*[0,ti] is mutually 
absolutely continuous with respect to the law of B^^ , with density proportional 
to l/Li{B) relative to the law of B^'^ . That is, for all non-negative measurable 
path functionals F 

where c is determined by the normalization of local time via Hll|) . 

While the laws of iJjf' [0, Du] and the pseudo- bridge i?*[0, n] are mutually 
singular, their random occupation measures have the same distributions. In 
fact, the sample path of B^^lQ, Dij] is simply a random rearrangement of the 
sample path of i?, [0, ti]: 

Lemma 13. Let U be a uniform (0, 1) variable independent of B^' , and 
independent of X distributed like B^[0,ti]. Then a process Y distributed 
like B^'^lO, Djj] is created by the following rearrangement of the path of X, 
whereby the random occupation measures of X and Y are pathwise identical: 
let {Gu,Dij) be the excursion interval of X straddling time U, and let Y be 
derived from X by swapping the order of the path fragments X[Gu, Du] and 
X[DuA], say 

Y ^X[0,Gu]:X[DuA]--X[Gu,Du] (31) 
with an obvious notation for concatenation of path fragments. 

Proof. By construction, the path of Y ends with a _B-excursion of length 
1 - Gi{Y) = Du - Gu- The joint law of y[0,G'i(r)] and Y[Gi{Y),l] := 
X[Gu^Du] was described in [211 Theorem 1.3] and 25, Theorem 3.1 and 
(3.d)], and is identical to the joint law of Z[0, Gi(Z)] and Z[Gi{Z),l] for 
Z B^'^[0, DuiB^")], which can be read from LemmaEl To be explicit, the 
common distribution of r[0, Gi{Y)] and Z[0, Gi{Z)] is that of B[0, Gi{B)] de- 
scribed by (|aj, while both y,[Gi(y),l] ■.^X4Gu,Du] and Z,[Gi(r),l] := 
B^'^lGuiB^'^), Dij{B^'^)] are standard _B-excursions. Since the excursion is in 
each case independent of the preceding fragment, it follows that Y ~ Z. □ 

Proposition 14. Fix ^ > 0. Let G be a random variable independent of B^'^ , 
with G = Fa/^. The distributions of the two bivariate sequences, defined by 



F{B'='') 



18 David Aldous and Jim Pitman 

the lengths and bridge local time measures of intervals of the D-partition and 
the T -partition respectively, are determined as follows: 

(i) For Ij — Ij^ the bivariate sequence 

{GXi,,G-lY^)^^^^ (32) 

is the sequence of points {Xj,Yj), in X-biased random order, of a Poisson 
process on ]R?,q with intensity measure 

v{dt, dt) at-^e^^' dt Pif^T^" E dl) = £-^P{n e dt)e-^' (33) 

for Ti as in ()ll|l. which makes 

S,X,^^ and S,Y,^^. (34) 

(ii) // the points {Xj,Yj) of a Poisson process with intensity v on R?,q are 
listed in X-biased order then 

for Sx —Ej^j as m (El). 

(iii) For Ij = the bivariate sequence in (|32|) is the sequence of points, 
say {X'j,Yj), in Y' -biased random order, of another Poisson process on R?,q 
with the same intensity measure v. So if in (ii) the points {Xj,Yj) are listed 
instead in Y-biased order, then H35II holds with the sequence of T -intervals 
instead of the sequence of D-intervals. 

Proof. Part (i) is proved in Section^] Part (ii) is just a restatement of part (i). 
Part (iii) is proved in Section □ 

Note that the normaUzation in (|35(> involves Ux and its ath power, both 
for the /^-partition and for the T-partition. Obviously, this is easier to handle 
if the sampling is X-biased rather than F-biased, which is one explanation of 
why various distributions associated with {Ij') are simpler than their coun- 
terparts for {if). 



4 Analysis of the D-partition 

As a preliminary for the proof of Proposition ^] (i) , we recall the following 
well known lemma, which characterizes the distribution of a sequence {Qj), 
known as the GEM{9) distribution after Griffiths, Engen and McCloskey. 
The distribution of {Q(j)) obtained by ranking {Qj) is known as the Poisson- 
Dirichlet distribution with parameter 9. See jl7) . [181 §9.6], |32| . 
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Lemma 15 (characterizations of GEM (9) |53J EH EHI ) . Fix 6 > and 

^ > 0. Let G and Qj,j = 1, 2, . . . be non-negative random variables. Then the 
following are equivalent: 

(i) the sequence (Qj) admits the representation Qj = Wj Y[i=ii^~^i) where 
the Wj are independent heta{l,6) variables, and G is independent of (Qj) 
with G = Fe/S,; 

(ii) ^jQj — 1 a.s. and (GQj) is the sequence of points of a Poisson point 
process on ]R>o with intensity Ot^^e^^* dt, listed in size-biased order. 

The next well known result 18, §5.2], 37, Prop. 4.10.1], combined with the 
previous lemma, provides an efficient way to identify various Poisson processes. 

Lemma 16 (Poisson marking). Let {S,S) and {T,T) be two measurable 
spaces. Let (Xj) and (Yj) be two sequences of random variables, with values in 
S and T respectively, such that the counting process S ■) is Poisson 

with intensity measure fi on S , and the Yj are conditionally independent given 
(Xj), with 

PiY,e-\x^,x,,...) = p'{x,,-) 

for some Markov kernel P' from {S,S) to {T,T). Then the counting process 
'^j ^{{Xj,Yj) (z ■) is a Poisson process on the product space S xT with 
intensity measure fi{dx) P' {x , dy) on the product a-field. 

Proof of Proposition \14\ (i). Proposition 1101 (i) and Lemma |15I (i) show that 
(Xjn,j ^ 1) has GEAIia) distribution. By assumption, G is independent of 

this sequence with G = Pa/S,. Lemmall^limplies that {GXjd ) is the size-biased 

j 

ordering of a Poisson point process of intensity at~^e~^'^ dt. PropositionllOlfii') 
and Lemma [T51 now identify the {GXjd jG^L^o) as the points of a Poisson 

process with intensity measure v defined by the first expression in (|33|l . To 
check the equality of the two expressions in let 

fi{t) := P{n e dt)ldt. (36) 

Since n = ^^/"n by 

fi{i)=r^' ''flit (37) 

whereas by another change of variables 

P(rTf " e dt)/dl = a-^tr^-^''^fi(t/e^''') = a-hr^fe{t) (38) 

and the identity follows. By application of Hll|l . the ^-marginal of i/ is 
£^^e~^'» ^ d£. The distribution of J^j infinitely divisible law with 

this Levy measure, that is the exponential distribution with rate c^". □ 

Implicit in Lemma IT^ and H38|l is the following formula of |25[ (3.u)] for 
the density of L° := L\{B'°'') 
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P{L° e d£) = car{a)iP{T-" £ d£) = cr{a)fi{l) dl (39) 

for fe{x) as in H36|l the stable (a) density of determined by Hll|l . That 
is to say, the distribution of L\{B^'') is obtained by size-biasing the com- 
mon distribution of Tj~" and Li{B). In particular, the general formula (|39|l is 
consistent with Levy's well known formulae in the Brownian case j21j, with 
a = l/2,c= V2 

/,(a;) = ^a;-3/2e-5^'/^ (40) 

and 

p{Llede)/de^ee-^^\ (41) 

For general a, a series expression for fe{x) is known [311 031 El EDI • If a = ^/n 

for some n = 2, 3, . . ., integral expressions for fi{x) can be derived from a 
representation of 1 /ti as a product of n — 1 independent gamma variables [431 
Theorem 3.4.3]. To conclude this section, we record the following immediate 
consequence of Proposition ^1 (i): 

Corollary 17. Let (aj) be a sequence of i.i.d. copies of ti with the stable 
(a) distribution and let (Qj) with GEM{a) distribution of Lemma \15\ 

be independent of {(Jj). Let Lj := [Qj/cfj)" and L :— Lj. Then L = L\ 
as in ()39|l . and the sequence (Lj/L) has GEM{1) distribution, independently 
of L. 



5 Analysis of the T-partition 

We start by recalling the structure of a Markov process up to the last time it 
visits its initial state before an independent exponential time. This does not 
involve the self-similarity assumption. 

Lemma 18 (jl4j). Let (t£,£ 0) be a drift-free subordinator which is the in- 
verse of the continuous local time process {Lf{B), t ^ 0) of a regular recurrent 
point 0, for a strong Markov process B started at 0. Let e be an exponential 
variable with rate ^, with e independent of B, and let 

G:^Ge{B) and L -.^ L%{B) = LI{B). (42) 

(i) The local time L has exponential distribution with rate the Laplace 
exponent of the subordinator defined by E{e~^'^') = e~'^^^^^ . 

(ii) For £ > 0, there is the equality in distribution of path fragments 

(B[0, G]\L^e)^ (B[0, Ti] I Ti < e). (43) 

(iii) The joint distribution of (G, L) is 



P{G edt,Le di) = V(0 d^e~«*P(r^ e dt). 



(44) 
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which is the distribution of the value at time 1 of a drift free bivariate subor- 
dinator with Levy measure 

v{dt,d£) = e-^d£e-^*P{Te e dt) (45) 

whose i-marginal is the Levy measure £~ie~^^^''^ d£ of the exponential distri- 
bution of L. 

Proof. These results are derived from Ito's theory of excursions of B, by letting 
{Nt,t ^ 0) be a Poisson process with rate ^, independent of B, and taking 
£ to be the time of the first point of N. To briefly recall the argument, say 
that a jump interval {Ty^,Ty) of the inverse local time process r is marked 
if N{Ty-,Ty\ > and unmarked otherwise. Then, by basic theory of Poisson 
point processes, the sum of unmarked jumps 

^" ■■= E - ry-)HN{Ty.,Ty] = 0) (46) 
o<y<e 

defines a subordinator with distribution 

F(t^" e dt) = c^'^^^^-^'Pin e dt). (47) 

The left end G of the first marked interval is G = = t^, and the subordi- 
nator summing unmarked jumps of r is independent of L, the local time 
of the first marked jump. See also [H EH ESI ■ □ 

To be more explicit, part (iii) of the Lemma states that 

iG,L) = iS^X^,EjYj) (48) 

for {Xj,Yj) the points of a Poisson point process on M?,q with intensity mea- 
sure v defined by (|45|) . In particular, for a self-similar B as in Section El this 
measure v is identical to the measure i/ featured in H33I) . 

In the setting of Lemma [THl even with construction of the Poisson process 
of marks of rate ^ independent of B, more randomization is required to con- 
struct points {Xj , Yj ) such that (|48|1 holds with equality almost surely rather 
than just in distribution. But this can be done by the following construction, 
which is the basis of our proofs of Proposition ^] and Proposition 1141 fiii'l . 

Lemma 19. In the setting of Lemma \lSl let I :— [Gj, Dj] be a random subin- 
terval of [0, 1), where the endpoint 1 is deliberately excluded, to avoid the 
jump of the inverse local time process (ti) at time £ — L in the following con- 
struction. Suppose I is independent of B and e, and define further random 
intervals 

IL := [GiL,DiL] and t{IL) [t{GiL),t{DiL)] (49) 
where t{£) := for £ ^ 0. 
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(i) For y>0, there is the equality in distribution of path fragments 

{B[t{IL)] I \jL = y) ' {B[0, Ty] \Ty<e) (50) 

where Xil '■= DjL — GiL is the increment of local time of B over the time 
interval t[IL). 

(ii) // {Ij) is an interval partition of [0,1) which is independent of e and 
B, and {tj-l) is the corresponding interval partition of [0,G), then given 
the sequence of local time increments (Xijl) the path fragments B[tj.]^] are 
conditionally independent with distributions described by (|50|l for Ij instead 

(iii) // Ij :~ [Vj-i,Vj] with Vj := 1 — ni=i(l ~ ^i) f^''' independent 
uniform(0, 1) variables Ui independent of B and e, then the bivariate se- 
quence of local time increments and path fragments 

{Xj^l,B[t{I,,L)])^^^^^_ (51) 

is the sequence of points of a Poisson point process on M>o x f2, in local-time 
biased order, for a suitable space of path fragments f2 of arbitrary finite length, 
with intensity measure 

£~^PiTt, < e, B[Q, Tf] e duo) (52) 

whose £-marginal is the Levy measure £~^e~''^^^^^ d£ of the exponential distri- 
bution of L with rate ipiO- 

(iv) Let Xj :— t{Gi-L) — t{Gi.L) be the length and Yj :~ XijL the local 
time increment associated with the random subinterval T{IjL) of [0, G). Then 
the (Xj,Yj) are the points of a Poisson point process on IR?,g with intensity 
measure v defined by (|45l) . in Y -biased random order, and 

{G,L) = [SjXj, SjYj) almost surely. (53) 

Proof. The first two assertions are straightforward consequences of the pre- 
vious Lemma. Part (iii) foUows from (ii), the Poisson representation of 
GEM{1) in Lemma [T5I and Poisson marking. Part (iv) foUows from (iii) and 
Lemma ^1 (iii) . □ 

Proof of Proposition^^ (iii). We will exploit the foUowing construction of the 
standard bridge B^^ of the self-similar Markov process B by random scaling, 
as in PI] and 35, Lemma 4]. Let 

B'^'' := B^ [0, Ge{B)] for e independent of B with e = AC, (54) 

so £ is exponential with rate S^. Then 

G:=G,{B)^^ and L := L%{B) = G'^ LI{B^') = ^ (55) 
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by H21|) . H19|) . and a-scaling of local times, where the exponential distribution 
of L is read from Lemma |18I Suppose now that Ij := [Vj^i,Vj], for Vj as in 
Lemma [T^ The T-sequence is now constructed as a function of these Vj and 
B^'^ := B^[0,G] as in according to Definition [21 that is 

Tj := M{u : L"JLBi >Vj}. (56) 

ByEHland 

Xi,L = m - = {V, - Vj-i) G'^LliB'-') - G'^L'>jj{B'"^) 

Part (iii) of Proposition E| can now be read from Lemma [T^ fivV □ 

Proof of Proposition^^ Parts (i) and (ii) follow immediately from the result 
of Proposition 1141 (iii) proved above. Turning to consideration of the path 
fragments, we observe by switching identity (|17|l that 

(i?, [0, n] I r, = = {B""' I L° = It"^) (57) 

where L\ := LKB^') as usual, and a regular conditional distribution for B^'^ 
given Li can be as constructed in :28:. Lemma 12]. Hence from (|52f) . if /?i 
denotes a suitable space of paths of length 1 , the trivariate sequence of local 
time increments, lengths of path fragments, and standardized path fragments 

{Xi,L,Xr,^,,B,[Ti^L])^^,^,_ = {G'^L°j{B'^'),GXjT,B'^^[lJ])^^^^^^ (58) 

is a Poisson process on R>o x M>o x whose intensity measure is 

d£P{n e dt)e~^'P{B^'' e dui \ L\ = ir"'). (59) 

Using the first form of v in H33|l to integrate out i in H59|l . we see that the 
lengths and standardized fragments 

(H.'^*[-/,i]),=i,2,... = (GAzj,i?!?^[/J]),.=i,2,... (60) 
form a Poisson process on R>o x Qi whose intensity measure is 

ar^Q-^^ dtQ{AuJi) (61) 

where 

/■oc 

g(dcc;i) = / F(B'^'' e dc^i I = y)P(Tf " e dy) P{B, [0, n] e diu) (62) 

by the switching identity (|57|l . The factorization in H61|l shows that the B^^[lJ] 
are i.i.d. copies of _B*[0, ri] when listed in length-ranked order. That is part (iii) 
of Proposition ^] □ 
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5.1 Further distributional results 

We record in this section a number of further formulae related to the distri- 
bution of the lengths and local times defined by the T-partition. 

Proposition 20. For the T-partition the (2n + l)-variate joint density of the 
total bridge local time L\, the lengths of the first n intervals, and the local 
times at on these intervals, is given by the formula 

p{lI e d£, \q e Axj^L'^'iiJ) e dy,, 1 s; i ==; n) 

= cr{a) difi-y^ -xi Xn) [[ 



j=i^-y^ yj-i 

for fy(x) := P(Ty G da;)/da; the stableia) density as in ()36|l . 

Proof. This follows from the switching identity H17(l and the definition of the 
T-partition. □ 

While the distributions of the cut times Tfc and interval lengths (Tk — Tk-i) 
in principle determined Proposition 1201 formulae for these distributions are 
more easily obtained as follows. For < u < 1, let 

T^^ M{t : LO/L? = u} 

where (L^ , ^ t ^ 1) is the local time process at of B^^'. Then by use of the 
switching identity ((T7|l we can write down for < a; < 1, < ^ < oo, 

P(r„b'- edx\L'i = £) /dx = /"^(^)/'«^(^) (63) 

/HI) 

where x := 1 — x. Integrating out with respect to the distribution of 
gives the density 



P(r^' e dx)/da; = cr{a) 

Jo 



fue{x)fue{^)de. (64) 



which can be simplified using Levy's formula H4U|) in the Brownian case to 
give for a = 1/2 

P(t„^'- e dx)/dx = 2T -2^3/2 (0 < < 1). (65) 

In particular, for u = 1/2 we recover the the result of 9, Theorem 3.2] that 
T]y2 has uniform distribution on [0, 1] in the Brownian case. 

According to Definition |21 Tk r-/, for Vk independent of B^'^ with 
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for iJfc is a product of k independent uniform(0, 1) variables, with 

mt^^tpi^ and V^(^Zi£M.i (66) 

because logTTfc is the fe'th point of a rate 1 Poisson process on [0, cx)). Since 
the process (t^J , ^ u ^ 1) has exchangeable increments, we find that 1 — Tk 
and the length of the fcth T-interval have the common distribution 



P(l - Tfc e dx) = P{XjT edx) ^ P{t^' e dx)P(i7fc G di 



(67) 



In particular, in the Brownian case a ~ 1/2, H67|) and (|65|l yield the curious 
formula 

P{Tiedx) h(x) + h(x) 1 , / I A 

Corollary 21. The point process of lengths of T -intervals has mean density 
Vp(Ajt edx) =aa;-i(l-a;)"-Ma;= / P(t^'' G da;)^- Mm (69) 



fc=i 



/or x e (0, 1). 

Proof. The first equality is read from part (i) of Proposition II II and the well 
known formula for the mean density of points of a Poisson-Dirichlet(Q;) dis- 
tributed sequence |82l (6)], which can be read from Lemma IT^ The second 
equality is then read from l|H7| and □ 

For general a, the second equality in (|69|) does not seem very obvious 
from H63I) and (|64|l . However, it can be checked for a ~ 1/2 using H65|) . and 
it can also be verified by a very general argument, which we indicate in Sec- 
tion 

Path decompositions of B^'^ at the times Tk are more complicated than the 
corresponding decompositions for the times Dy. expressed by Lemma |H1 For 
the T-partition, the pieces are not pure _B-bridges. Rather, when normalized 
they have density factors involving their local times at 0. Compare with similar 
constructions in [1111131 I^l33j . 

By the Poisson analysis of the previous section, conditionally given (Ti , 
LSJn^, LJ) the pieces of P^'' before and after time Ti are independent P-bridges 
with prescribed lengths and local times at 0. The appearance of ft. -I- fc in 
formula (a) below shows that the right side does not factor into a function of 
{x, h) and a function of {x, k). So even in the Brownian case, L^^ and L\ — L^_^ 
are not conditionally independent given Ti, and hence the same can be said 
of the fragments of B^^ before and after time Ti. 
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Proposition 22. In the Brownian case with a = 1/2, c = \/2, 

F(ri e da;, e d/i, L? - e dfc) hk (a;x)-3/2 / P 

— ■ exp ' 



dxdhdk V27r {h + k) \ 2a; x 

while for 

X:^^= L\{Bl^%m and Y ^^ZJ^ = L\{B'1^{T^, 1]). 

i/iere is t/ie joint density 

P(Xeda,red6) a6 / a" 

/(a,5)exp -— -— (70) 



dad5 V27r V 2 2 

w/iere 

/(a b) /' J^^ll^dx^ilogf^^^i^^^^ 

io a^i: + r \{r-a){r-b)^ 

for r := va? + &^. 

Proof. The first formula is an instance of Proposition l20l which we now check. 
With notation as in (|^ . 

P(ri e da; I L"j,^ =h,Ll- 1°^, = fc)/da; = f{x \ u, £) 

for h — ui and k — u£. Wc also know, by definition of Ti, that 

P{L'^^ e d/i, - i?,^ e dfc) = d/i dfc e" 

where ^ = h+k and an has canceled the factor of £ in the density (|41|l of i". 
Combining these formulae gives the trivariate density of (Ti , Lj,^ , Lj* — L j.^ ) , 
which rescales to give 

P(Ti £ da;, X g da, Y G db) ab jxx)-^/^ f ^ ^ 

da; dad6 ~ {a^ + by/^) ^^"^^Y ^ T 

and (|70() follows by integrating out x. □ 



6 Complements 

6.1 Mappings conditioned to have a single basin 

In the Brownian variation of the transformation from X to F in 

Lemma which further swaps the exchangeable pair of fragments X[Q, Gu] 
and X[Du, 1], is the continuous analog of the transformation, mentioned in 
Fact H2.2|l fe) from the stretch of the cycles-first mapping walk for a given 
basin to the stretch of the basins-first walk for the same basin. As pointed 
out in the last section of pP, if the uniform mapping of [n] is conditioned to 
have only one cycle, the scaled basins-first walk converges in distribution to 
the process 2\B]^'^[0, Di/{B^'^)]\. The above argument yields: 
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Corollary 23. For a uniform mapping of [n] conditioned to have only one 
cycle, the scaled cycles-first walk converges in distribution to 2|_B»[0, ti]\ where 
i?*[0,Ti] is the Brownian pseudo-bridge. 

The distributions of several basic functional of pseudo-bridge [0, ri] 
are known. In particular, the occupation density of the reflected process is 
governed by the same stochastic differential equation governing the occupation 
density process of a reflecting Brownian bridge or Brownian excursion |28j . 
According to Knight ^H] (see also ^^'^ papers cited there) , the law of the 
maximum of the reflected pseudo-bridge is identical to that of l/(2-\/_ffi(i?3)) 
where Hi^R^) is the hitting time of 1 by the three-dimensional Bessel process, 
with transform E{c^p{-^e'^Hi{R3))) = 6'/sinh6' for real 9. Thus we deduce: 

Corollary 24. For a uniform mapping of [n] conditioned to have only one 
cycle, the asymptotic distribution of the maximum height of any tree above the 
cycle, normalized by \fn, is the distribution of \/ Hi{R^). 

See also ^21 for a survey of closely related distributions and their applica- 
tions. 

6.2 Exchangeable interval partitions 

Suppose that {If^) is an exchangeable interval partition of [0, 1]. That is (as- 
suming for simplicity that the lengths A/ox are almost surely all distinct), for 
each n — 2, 3, . . . such that A/ox^ > 0, where (-^(°j)) is the associated length- 
ranked interval partition, conditionally given A/ox^ > the ordering of the 
longest n sub-intervals 1°^-^ , 1 ^ j ^ n is equally likely to be any one of the n\ 
possible orders, independently of the lengths of these n intervals. Call {I'j^) 
infinite if ^(A/ox^ > 0) = 1 for all n. As shown by Kallenberg jJS], for an 
infinite exchangeable interval partition {If^)^ for each u £ [0,1] the fraction 
of the longest n intervals that lie to the left of u has an almost sure limit L° 
as n — > oo. The process (i",0 ^ u ^ 1) is a continuous increasing process, 
the normalized local time process of {Ij^)- It is easily shown that for B^^ 
as in previous sections, and more generally for B^'^ the standard bridge of 
any nice recurrent Markov process, constructed as in |13| . the interval par- 
tition (Ij^) defined by the excursions of B^'^ away from is an infinite ex- 
changeable interval partition of [0,1], whose normalized local time process 
is — L'^/LijO ^ u ^ 1 for any of the usual Markovian definitions of a 
bridge local time process := L^(_B'"'). In particular, this remark applies to 
a self-similar recurrent process B as considered in previous sections. 

Theorem 25. The assertions of Theorem\^ remain valid for the D- and T- 
partitions defined by Definitions^ and\^ for any infinite exchangeable inter- 
val partition {Ij^) instead of the excursion intervals of a standard Brownian 
bridge B^^ , with the complement of [J,- 1°^ in [0, 1] instead of the zero set of 
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B^^ , and the normalized local time process (L" , ^ u ^ 1) of (Ij^) instead of 
{L'^/Li,0 ^ u ^ 1). Moreover, the sequence of normalized local times {L^i.) 
has the same GEM{1) distribution for Ij = Ij^ as for Ij — ij . 

Theorem|2nican be derived from a certain combinatorial analog, stated and 
proved as Lemma 1^ below. Let us briefly outline the method of derivation, 
without details. Consider an infinite exchangeable interval partition {Ij). Take 
k independent uniform (0, 1) sample points, assign "weight" 1/fc to each, and 
let {if^^) be the intervals containing at least one sample point. Each interval 

I^^^ is thereby assigned weight 1/fcx (number of sample points in interval). 
For fixed k we can apply Lemma interpreting "length" as "weight", and 
conditionally on the number of intervals in the partition. The conclusion of 
Lemma |26l is a variant of the desired Theorem 13 for (/j ), in which "position 
X G (0, 1)" of interval endpoint is replaced by "1/fcx (number of sample points 
in (0, x)y\ and in which "normalized local time at -u S (0, 1)" is replaced by 
"relative number of sampled intervals in (0,m)". One can now argue that as 
/c — > cx) we have a.s. convergence of these variant quantities to the original 
quantities in Theorem 1251 

Lemma 26. Let (/°'^)i^i^n be an exchangeable interval partition of [0, 1] into 
n subintervals of strictly positive length. Define Dy as in Definition Q for 
1 ^ J ^ , where Jn is the first j such that Dy ^ 1, to create a D- 
partition {l[')i^j^jo of [0,1], and define a T-partition {lj)i^j^j'r of [0,1] 
similarly using cut points Tj,l ^ j ^ determined as follows: given that the 
random set of endpoints of {Ij^)it^js^n is lA := {ujjo^j^n with = uq < ui < 
■ ■ ■ < Un — 1, let Ti have uniform distribution on hi D (0, 1] , and given also 
Ti — ti < 1 let T2 have uniform distribution on U D {ti, 1], and so on, until 
TjT — 1. For Ij an interval of either of the D- or T -partitions so defined, let 
Ni- denote the number of intervals of (/f )i^i^ri which are contained in Ij, 
so 1 ^ Nj. ^ n. Then the assertions of Theorem\^ remain valid provided that 
Nj^ is substituted everywhere for iPj, . 

Proof. We will check that part (i) of Theorem holds in this setup, along 
with 1)72(1 . The remaining claims are straightforward and left to the reader. By 
conditioning on the ranked lengths A/jjx of the intervals (/f'^)i^i^n, it suffices 
to consider the case when these ranked lengths are distinct constants. Let 
denote the random partition of [n] defined by the random equivalence 
relation i '--^ j iff and I°J^ are part of the same component interval of 
the Z)-partition, and define TTj similarly in terms of the T-partition. Since 
each unordered collection of lengths and sub-interval counts is a function 

of the corresponding partition, it suffices to show that 77^ = TTJ. Due to 
the well known connection between the discrete stick-breaking scheme used 
to define the T-partition and the cycle structure of random permutations, 
which was recalled in Section (d), we can write down the distribution of 
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TTJ without calculation: for each unordered partition of [n] into k non-empty 
subsets {Ai, . . . , Afc}, 

k 

P{n^ = {A„..., A,}) = 1 J] (1^^. I - 1) ! (71) 

^' i=i 

where \Ai\ is the number of elements of Ai. On the other hand, for the D- 
partition, for each ordered partition (Ai, . . . , Ak) and each choice of aj G 
Aj,l ^ j ^ fc, with A(a) the length of I^^^ and X{A) :— J2aeAM'^)j can 
write down the probability 

P(lf-[J /(X) has right end Dy^ ^ ^11(1^,1-1)! ^/^"j, . 



where the factors of — 1)! come from the different possible orderings of all 
but the last I^^ to form ij^ . If we now sum over all choices of aj S Aj , for each 
i ^ j ^ k, we find that X{aj) is simply replaced by X{Aj). If we then replace 
(Ai, . . . , Aj.) by (Act(i), . . . ,v4.(o.(fc)) and sum over all permutations a of [fc], 
to consider all sequences of sets consistent with a given unordered partition 
{Ai, . . . , Afc}, we get precisely (|71ll for 77,f instead of TTj, due to the identity 



k 



En 



^k , ^ 



= 1 T,t=j H^'rit)) 

This is obvious, because the product is the probability of picking the sequence 
of sets (Ao.(j), I ^ j ^ k) in a process of A(Ai)-biased sampling of blocks of 
the partition {Ai, . . . ,Ak}. □ 

We note the consequence of the previous proof that the number of com- 
ponents of the I?-partition and the number of components Jj of the 
T-partition have the same distribution, which is the same for every exchange- 
able interval partition {If^)i<^i<^n of [0, 1] into n subintervals of strictly positive 
length: 

n 

I Jj ^K,,^Y: ^cu (72) 

4=1 

where Kn is the number of cycles of a uniformly distributed random permuta- 
tion of [n], and the Ci are independent events with P{Ci) = The second 
two of these equalities in distribution are well known and easily explained 
without calculation "26". But the first is quite surprising, and we do not see 
how to explain it any more simply than by the previous proof. 



6.3 Intensity measures 



In this section we check Corollary 1211 bv showing it can be generalized and 
proved as follows: 
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Corollary 27. In the setting of Theorem \25l the common intensity measure 
of the the point process of lengths of T -intervals and the the point process of 
lengths of D -intervals is 

k=l fc=l 

(73) 

where '■= iid{t : > u} is the inverse of the normalized local time process 
of the exchangeable interval partition. 

Proof. The equahty of the first three measures displayed in (|73|l is read from 
the conclusion of Theorem l25l using the fact that the ZJ-partition is in length- 
biased order. The equality of the first and fourth measures follows from the 
definition of the Tfc, the exchangeable increments of (?„, ^ u ^ 1), and (|66|l . 
just as in the proof of 1)69(1 . □ 

As a check on Theorem 1251 let us verify the equality of the second and 
fourth measures in ((73|l in the following special case, which includes the setting 
of Corollary ED 

Let {ti,£ ^ 0) be the inverse local time process of B at 0, for B as in 
Lemma [THl not necessarily self-similar. Note that we must explicitly assume 
{t(,£ ^ 0) is drift free for the conclusion of part (iii) of that Lemma to be 
true. We assume that now. Assume that the Levy measure of (r^,^ ^ 0) has 
density p{x). Let {Ij^) be the exchangeable partition of [0, 1] generated by the 
excursion intervals of B conditional on i?i = and Li{B) — I for some fixed 
£ > 0, or equivalently by the jumps of (ts,0 ^ s ^ £) given ti = \. Then, 
formula H63|l generalizes easily to show that the fourth measure in H73|l has 
density at x 



(74) 

for fi{x) as in H36|) . On the other hand, abbreviating D :— Dv^ and G :— Gvi 
so [G, D] is the interval I'^^ which covers the independent uniform time Vi , 
we know from (|75l) that for < w < 1 

P(l -{D-G)^ Aw) = ^P^' - - ^)^Md^ . 

Also, it is easily seen that conditionally given 1— (D — G) = w, the normalized 
local time Lq is uniform on (0, 1) and independent of the pair (G, 1 — D), which 
is distributed like (t„^, Tui) conditioned on Te = w. Together with the previous 
formula for w = y -\- 1 — x^ this gives the trivariate density 

P{L''(.edu,Gedy,Dedx) ^ ep{x-y){x~y)fuiiy)fuiix) 
dudydx fe{l) 

{0 < y < X <1). Now lO implies that 
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fui{y)p{x -y)[x- y) dy 



ui 



so we deduce that 




and hence that the density displayed in (|74|l is indeed x ^P{D £ da;)/da;. 
6.4 Two orderings of a bivariate Poisson process 

According to Proposition^] for each a G (0,1) the Poisson point process with 
intensity measure v(dx, dy) = p{x, y) dx dy displayed in (|33|l has the following 
paradoxical property: 

(a) If the points {Xj^Yj) are put in X-biased order, then the Yj are in Y- 
biased order, whereas 

(b) if the points {Xj^Yj) are put in F-biased order, then the Xj are not in 
X-biased random order; even the distribution of Xi is wrong. 

We first see this for a = 1/2 by passage to the limit of elementary combinato- 
rial properties of uniform random mappings. We then see it for general a from 
the bridge representations of Proposition^] Other point processes of lengths 
and local times with these properties can be constructed from an exchange- 
able interval partition, as shown by Theorem 1251 in the previous section and 
Lemma [T51 This argument, shows that (a) holds for the bivariate Poisson pro- 
cess with intensity (|45|l featured in Lemma IT^ for any drift free subordinator 
(Ty,y ^ 0) with E(e~^^^) = e""^*^^-*^. Then the Yj normalized by their sum 
have GEM{1) distribution, both for an X-biased and for a y-biased ordering. 
We offer here a slightly different explanation of (a) in this case. That is, given 
some joint density p{x,y), we indicate conditions on p which are necessary 
and sufficient for (a) to hold for the bivariate Poisson process with intensity 
p, and then check that these conditions are in fact satisfied in the case (|45f) . 

Let {Xj,Yj) be the points of a Poisson process on R^q with intensity 
p{x, y) dx dy, in X-biased order. Let Ex '■= X]j -^j ^^'^ '■= Yj- Let 



fx{x) -.^PiSx e dx)/dx; 



fviy) 



P{Uy e dy)/dy 




By a basic Palm calculation, as in 



P{Xi e da;, Sx "Xie dw) 



px (x) dx fx (w) dw 



X 



(75) 



X + w 



and similarly, with 
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fx,Yix,y) := PiSx e dx,SY e dy)/{dxdy) 



P{Xi e da;, Yi G dy, Sx - Xi e dw, Sy - Yi e dv) 



p{x, y) dx dy fx.Y (w, v) dw dv 



X + w 



(76) 



Now, a necessary condition for the Yj to be in ^-biased order is that Yi should 
have the same joint distribution with Uy as if Yi were a size-biased pick from 
the Yi, that is like {TSJ 



P(Yi e dy, i:y - Fi e dv) - (y) dy /y (w) dv 



y 



y + v 



(77) 



Thus a necessary condition on p{x, y) for (a) to hold is that for all y, u ^ 



dx 



dw p{x,y)fx.Y{w,v) 



X 



X + w 



PY{y)fY{v) 



y 



y + v 



(78) 



Moreover, by keeping track of the first k of the {Xj, Yj) jointly with 2Jx and 
Sy it is clear that we can write down a multivariate version of (|78|) whose 
truth for all k would be necessary and sufficient for (a). 

In the special case (1^ . with fy{x) :— P{Ty £ da;)/da;, the subordination 
argument gives 

pix,y) = y-^ fy{x)e-'^'' . 
Since the Y-marginal is exponential with rate 

fY{y)^mypY{y) = me-''^^^'' 

and hence by generalization of (|43|l . using {Sx, ^y) = (G, L), 

fx^Y{x,y) = ^{Ofy{x)e-^''' = ■ip{^)yp{x,y). 

If these expressions are substituted in l|78|) . and we use the definition of V'(C) 
on the right side, we find that (|78|l reduces to the identity 



E 



Ty + Ti,e-?('^!'+'^") 



y 



y + v 



But this is true by virtue of 
E 



y 



y + v 



which holds by exchangeability of increments of {ti,£ ^ 0). Moreover, the 
multivariate form of (j78(l mentioned above is easily checked the same way. 
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